arXiv: 1508.07073vl [math.OC] 28 Aug 2015 


1 


Minimum Number of Probes for Brain Dynamics Observability 
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Abstract 

In this paper, we address the problem of placing sensor probes in the brain such that the system dynamics’ are 
generically observable. The system dynamics whose states can encode for instance the fire-rating of the neurons or 
their ensemble following a neural-topological (structural) approach, and the sensors are assumed to be dedicated, 
i.e., can only measure a state at each time. Even though the mathematical description of brain dynamics is (yet) 
to be discovered, we build on its observed fractal characteristics and assume that the model of the brain activity 
satisfies fractional-order dynamics. 

Although the sensor placement explored in this paper is particularly considering the observability of brain 
dynamics, the proposed methodology applies to any fractional-order linear system. Thus, the main contribution of 
this paper is to show how to place the minimum number of dedicated sensors, i.e., sensors measuring only a state 
variable, to ensure generic observability in discrete-time fractional-order systems for a specified finite interval of 
time. Finally, an illustrative example of the main results is provided using electroencephalogram (EEG) data. 

I. INTRODUCTION 

In the recent years, control theoreticians have resorted to the tools available to study, analyze and design dynamical 
systems to better understand the dynamical behavior of the brain |[T1, |[2l, Q, Q, O, ||6l, O. Often, the dynamic 
model of the brain assumes a behavior at the neuron level, that is then analyzed in interaction with different neurons 
to which it communicates with lO, Q. A commonly used model for the neuron is the Hodgkin and Huxley, but 
several alternatives have been proposed. Most of the models have the following similarity: a neuron(i.e., beams 
a signal to neighboring neurons) if a certain threshold (or, gain) on its state (or stimulus) is reached. Nonetheless, 
the evolution of the state is believed to be non-linear in general, with a time-varying threshold. Therefore, several 
experimental efforts concluded that several regions in the brain exhibit specific properties. In particular, the non¬ 
linear dynamics in certain regions can be sequentially approximated by linear time-invariant systems, leading to 
a linear time-varying systems that can then be restricted to a finite possible dynamics due to some plasticity and 
invariant properties of the brain regions. Alternatively, the response to a stimulus may exhibit a long-range memory 
behavior which may indicate the presence of a fractional-order dynamics; for instance, the vesicular-ocular system. 
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the fly motion sensitive neuron and neocortical neurons |8j|, just to name a few. Nonetheless, several other dynamical 
systems are descrihed hy similar dynamics and sensing capahilities. For instance, fractional-order dynamics include 
artificial pancreas ||9l, on-chip traffic regulation lITOl . traffic dynamics across road arteries ifTTI . electrical and thermal 
machines, transmission and acoustics m, just to name a few. 

Hereafter, instead of determining a specific model to the brain’s dynamics given the neurons activity (or an 
ensemble of these) and their interactions, to which several studies have been conducted lIT^ . llT4l . we assume that the 
brain dynamics can be described by fractional-order dynamics. Then, we aim to provide a methodology to determine 
the minimum placement of dedicated sensors, i.e., sensors that measure a single state of the system, ensuring that the 
system is generically observable', more precisely, by generically observable (also known as structurally observable 


- see Section III for formal definition), we mean that almost all parametrizations (that are unknown) for a given 
structure of interaction (that can be known or imposed by neural-topological analysis). The reason we focus in 
determining such placement in the context of brain observability dynamics is threefold: (i) some probes (sensors) 
are invasive, i.e., require surgery to deploy the probes in the brain, which may be life threatening; (ii) nowadays 
deployment of probes in the brain neglects entirely the brain dynamics; and (iii) data collected in these probes is 
commonly used to feedback into the system, therefore, a smaller collection of data leads to a faster response time 
with potential implications in neurodegenerative diseases such as Parkinson’s and epilepsy ifTSl . 

To the best of the authors knowledge there are no systematic methods to do probe (sensor) placement in the 
brain. Instead, researchers have focused on collecting data from different regions of the brain and/or developing 
new sensors leveraging some of the new materials (some biodegradable that do not require the extraction after a 
period of time). On the other hand, if we consider dynamical systems, most sensor placement methodologies are 
either heuristic or greedy due to the combinatorial nature of the problem, see ifT^ and references therein. Structural 
systems theory studies the coupling between state variables rather than the specific strength of the coupling, suitable 
to the problems explored hereafter, since these are not available neither measurable directly, therefore they need to be 
estimated after data has been collected. The structural treatment of dynamical properties of switching systems and, 
in particular, the obtention of sufficient and necessary conditions to ensure structural observability, i.e., structural 
observable for some non-empty window of time, have been explored in ifTTl , iTT^ . This approach contrasts with the 
one presented in |[T^ . where conditions to ensure structural observability for all non-empty windows of time. In lIT^ . 
the deployment of dedicated sensors was explored, but the sensor placement to ensure structural observability for 
some non-empty window of time, has not yet been addressed. The present work also extends lIThl . Il20ll . where the 
minimum placement of dedicated sensors to ensure structural observability was addressed for linear time-invariant 
systems. o 

Main Contribution: The main contributions of the paper is to show how to place the minimum number of 
dedicated sensors to ensure structural observability in discrete-time fractional-order systems for a specified finite 
interval of time. 

The rest of the this paper is organized as follows. First, in Section |I^ we present the formal statements addressed 
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in the present paper. We introduce and revise some preliminary concepts and results in Section III In Section IV 


we present the main technical results, followed hy an illustrative example in Section |V] Finally, conclusions and 


discussion avenues for further research are presented in Section VI 


II. Problem Statement 

In this section, we formally introduce the dedicated sensor placement for discrete-time fractional order dynamics. 
Formally, this problem is described as follows. 

Vi Consider a model dynamics described by a linear discrete-time fractional-order system as follows 


A“xfc+i 




fc+i 

= '^AjXk, 
j=0 


( 1 ) 


A"" 


X 


n 

k-\-l _ 


where a* E R+{i = 1,... ,n), Aq = A and Aj = diag(-(-iy+i(^“\),-(-lV+i(^“\),...,-(-lV+i(^‘^"J), for 
brevity, we refer to the dynamics in ([T]) as F{A]ol,K) which is associated with the representation of all matrices 
Aj, with j = 0,..., K. Please note that a is a non-integer number and the term (“) = r(i+K-i+i) expressed 
via the Gamma function r(x) = 1121- 

Hereafter, given F{A\cx.,K) we aim to determine the minimum number of dedicated sensors J required to 
ensure that the linear discrete-time fractional-order system ([T]) with measured output 


yk = In^k, 


( 2 ) 


is structurally observable, i.e., 

arg min 1^71 

s.t. {F{A]cy., K),l^) is structurally observable, 

where J A {1,... ,n} denotes the set of indices of the dedicated sensors, represents the subset of columns of 
the identity matrix with indices in J, and {F{A-,a,K),l^) describes the system ([T])-Q. 

Notice that given K,K' where K' > K, we may have two different solutions to ([T]) with F{A]q.,K) and 
F{A]ol,K') for an arbitrary k. 


III. PRELIMINARIES AND TERMINOLOGY 

In this section, we review some notions of observability to discrete-time fractional-order systems, and their 
counterpart using structural systems theory ll22]| . 

We start by recalling that a solution in closed-form to ([T]) can be determined and given as follows. 

Lemma 1 ( /[27I/ . yf72l/ ).- The solution to ([T]) is given as follows: 


— Gk-\-lXo 


( 4 ) 
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where 


Gk 


A for /c = 0, 

^ ^^ 

x; AjGk-i-j for k>l. 

j=0 


o 


The notion of observability for discrete-time fractional-order systems is revisited in the next definition. 
Definition 1: The fractional-order system modeled hy ([T]l and Q is observable at time A: = 0 if and only if there 
exists some K > 0 such that the state xq at time k = 0 can he uniquely determined from the knowledge of the 
measured output yk with k = 1,..., K. o 

Of a particular interest is the notion of the observability matrix that is related with the system’s ohservahility. 
Definition 2 /TT?]/ ).' The ohservahility matrix associated with the discrete-time fractional order dynamics 

for a given time k as descrihed in ([T]), and the measured output 


Vk = Cxk, (5) 

is given as follows: 

ol = [iCGoy (GGiY ... {GGk-iYV. 

o 

Additionally, we have the following result. 

Theorem 1: The system descrihed hy ([T]l-Q is ohservahle if and only if there exists a finite time K such that 
rank (O;^) = n. o 

Given that the precise numerical values of the network parameters are generally not available for the large- 
scale systems of interest, a natural direction is to consider structured systems ll24l based reformulations of the 
above topology design problems, which we pursue in this paper. Representative work in structured systems theory 
may be found in ll25ll . lf2^ . |[27ll . lf28l . see also the survey and references therein. The main idea is to 
reformulate and study an equivalent class of systems for which system-theoretic properties are investigated based 
on the location of zeroes/non-zeroes of the state space representation matrices. Properties such as observability, in 
this framework, referred as structural observability, initially formalized for linear time-invariant systems ll25l . In 
fact, structural observability is a generic property, i.e., almost all (with respect to the Lebesgue measure) realizations 
satisfying a given structure are observable ll27l . Given the similarity between observability criteria of the linear¬ 
time invariant and the discrete-time fractional order dynamics, one can readily extend the notion of structural 
observability to the latter. More precisely, a pair {F{A;cx,K),G) is said to be structurally observable if there 
exists a pair {IF{A'; cx, K),G') with the same structure as {lF{A;cx,K),G), i.e., same locations of zeroes and 
non-zeroes, such that {lF{A'-,a,K),G') is observable. By density arguments ETl . it may be shown that if a pair 
{T{A;a,K),G) is structurally observable, then almost all (with respect to the Lebesgue measure) pairs with the 



5 


same structure as {F{A] q, K), C) are observable. In essence, structural observability is a property of the structure 
of the pair a, K), C) and not the specific numerical values. 

Now, we associate with each subsystem a directed graph (digraph) T> = T>{A, ) = (V, E), referred to as system 

digraph, with vertex set V and edge set E, where V = A Li y with A = {x^,... ,x'^} and 3 ^ = 
represents the state and output vertices, respectively. In addition, E = Ex,x U Exy where Ex,x = {{x^,x^) : 
^jk 7^ 0} and Ex,y = y-^) : j G J} represents the state edges and output edges, respectively. Similarly, 

we can define a state digraph T>{A) = {A,Ex,x)- Later, the matrix A is given in terms of A = Ai V ... V Am, 
where V corresponds to the entry-wise operation where if at least one of the entries is non-zero, then it provides a 
non-zero entry, and zero otherwise. A state vertex is said to be non-accessible by an output vertex if there exists 
no directed path, i.e., a sequence of directed edges where every edge ends in a vertex that is starting of another 
edges and no vertex is used twice, the state vertex to any output vertex. Additionally, we need to introduce the 
notion of a bipartite graph B{M) associated with a mi x m 2 matrix M given by B{M) = {Ti,C,Ec,n), where 
TZ = {ri,..., Tmi } and C = {ci,..., Cm 2 } correspond to the labeling row vertices and column vertices, respectively; 
further, Ec,ti = : Mij / 0}. The bipartite graph is an undirected graph with vertex set given by the union 

of the partition sets C and TZ, which we refer to as left and right vertex sets, respectively. A matching M C Ec,n 
is a collection of edges that have no vertices in common. A maximum matching is a matching with maximum 
cardinality among all possible matchings. For ease of reference, if a vertex in the left and right vertex set does 
not belong to an edge in a maximum matching we refer to as a right- and left-unmatched vertex, respectively. In 
addition, we can consider weights associated with the edges in a bipartite graph, so we can consider the problem 
of determining the maximum matching with the minimum sum of the weights, that we refer to as the minimum 
weight maximum matching. A digraph T>s = (V 5 , Es) is a subgraph of D = (V, E) if V 5 C V and Es C E. Finally, 
a strongly connected component (SCC) is a maximal subgraph (there is no other subgraph, containing it, with the 
same property) Vg = (V 5 ,£^ 5 ) of V such that for every u,v there exists a path from u to u and from v 

to u. We can create a directed acyclic graph (DAG) by visualizing each SCC as a virtual node, where there is 
a directed edge between vertices belonging to two SCCs if and only if there exists a directed edge connecting 
the corresponding SCCs in the digraph V = (V,f), the original digraph. The DAG associated with 'D{A) can be 
computed efficiently in 0{\V\ + |£’|) 1291 . The SCCs in the DAG may be further categorized as follows. 

Definition 3: lIThll An SCC is said to be linked if it has at least one incoming/outgoing edge from another SCC. 
In particular, an SCC is non-bottom linked if it has no outgoing edges from its vertices to the vertices of another 
SCC. o 

Finally, consider a mi x m 2 matrix M, and let M = {P G . p.^ = 0 if Mij = 0}, then the generic 

rank (g-rank) of M is given by g-rank(M) = max rank(P). 

PeM 

Now, we revise the structural observability necessary and sufficient conditions for linear switching systems. 

Theorem 2 f /f7^ ).- Consider a linear continuous-time switching system 
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x{t) = A^(t)x(t), (6) 

where a : M"'' M = {1,, m} is a switching signal, and x{t) G M” the state of the system at the instance of 
time t. In addition, let the measured output to he given hy 

y{t) = Cx(t). (7) 

The linear continuous-time switching system Q-Q is structurally ohservahle if and only if the following two 

conditions hold: 

(i) ^{Ai V ... V Am,^n) has no non-accessihle state vertex; 

(ii) g-rank{[Ai,..., Am,Ii[]) =n. o 

In addition, consider the following intermediate results. 

Lemma 2 ( HWS ): 'D{Ai V... V Am, ) has no non-accessihle state vertex if and only if there exits an edge to an 
output vertex in V{Ai V ... V Am,k^) from a state vertex in each non-hottom linked SCC of the DAG associated 
with T){Ai V ... V Am)- o 

Lemma 3 ( ^2^ ): There exists a maximum matching of B{M) with size n if and only if g-rank(M) = n. o 

IV. Main Results 

In this section, we present the main results of this paper. More precisely, we provide the solution to Vi that 
can he determined using Algorithm [T] Algorithm [T] is shown to he correct and with polynomial complexity in 
Theorem IH 

First, we provide necessary and sufficient conditions to ensure structural ohservahility of discrete time fractional- 
order dynamical systems. 

Theorem 3: A discrete time fractional-order dynamical system ([T])-Q is structurally ohservahle at time K if and 
only if the following two conditions hold: 

(i) T){Gq V Gi V ... V GK,"^n) has no non-accessihle state vertex; 

(ii) g-rank([G'o,Gi,...,Gif,I^]) = n, 

where Gi is as described in Q. o 

Proof: First, we notice that a linear continuous-time switching system Q-Q is structurally ohservahle if the 
ohservahility matrix associated with it given hy 

Oi = [G^ (CAi)T ... [GAmY (CAiA 2 )t (CAsA^t 

... ... a;’^))t .,, (cAri)T... {GA^fYYY 

has generic rank equal to n. From an algebraic point of view, if we set A^-^q = Gi-i (with z = 1,..., iT -|- 1) and 
G = lY■ Then, following the same steps as in Theorem 4 in lITSl . we obtain that the generic rank of equals 

that of O^. Consequently, conditions (i)-(ii) follow by invoking Theorem^ ■ 
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As a direct consequence of Theorem by invoking Lemma and Lemma we obtain the following result. 
Corollary 1: A fractional-order dynamical system ([T])-Q is structurally observable at time K if and only if the 
following two conditions hold: 

(i) there exits an edge to an output vertex in T>{Gq V Gi V ... V Gki^^u ) from a state vertex in each non-bottom 
linked SCC of the DAG associated with D(Go V Gi V ... V Gjc)> 

(ii) there exists a maximum matching of B{[Gq V Gi V ... V , I^]) with size n, 

where Gi is as described in Q. o 

ALGORITHM 1: Computing a minimal dedicated output placement ensuring structural observability 
Input: Gi, with i = I,... ,K 

Output: Minimal dedicated output placement 

Step 1. Determine the non-bottom linked SCCs Nf, i € I = {1, • • • , /3}, of T>{Gq \J .. .\J Gk) = (A, Sx,x) 

Step 2. Consider a weighted bipartite graph B{[Gq, ..., G]^, S']) = {TZ, C, Cnfi), where S is a n x /3 matrix and 

Sij = 1 if Xi G , and the column vertices be re-labeled as follows: the columns of Ai are indexed by {cl,..., c]j}, 
and the columns of S are indexed by {si,..., sp}. In addition, let the weight of the edges 

e S 7^ X y {cl,..., c)j} be equal to zero, the weight on the edges e GTZ x {si,..., s^} be equal to one, and 
all non-existing edges corresponding to inhnite weight. 

Step 3. Let M' be the maximum matching incurring in the minimum cost of the weighted bipartite graph presented in 
Step 2. 

Step 4. Take J' = {i (., s^) G M'}, i.e., the column vertices of S that belong to the edges in the MWMM M' (i.e., 

those with weight one). In addition, let J" = {1,...,n} \ {A: : (., Cj) G M',j = 1,..., n}, and 

J'" = U min{j :XjGAf^}. 

PG{1,...,/3}\J' 

Step 5. Set J = J' U J" U J"'. 


Now, we provide an efficient algorithmic procedure to compute a solution to Vi , described in Algorithm [T] that 
ensures that both conditions in Corollary [T] are satisfied. Briefly, Algorifhm [T] consisfs in finding a minimum weighf 
maximum mafching (MWMM) of a bipartite graph B{[Gl,..., G]^, S']), where the matrix S has as many columns 
as the number of non-bottom linked SCCs in the DAG representation of T>{Gq V... V Gk), and the non-zero entries 
in column i of S correspond to the indices of the state variables that belong to the f-th non-bottom linked SCC. In 
addition, we consider weights in the edges of the bipartite graph: those associated with the nonzero entries of Ai 
have zero weight, unitary weight is considered to the nonzero entries in S, and infinite weight otherwise. 

Therefore, if an edge with unitary weight belongs to the MWMM, then it contains a row vertex, which implies 
that dedicated sensors need to be assigned to the state variable with the index of the row vertex, whereas all the 
remaining unitary edges not in the MWMM mean that these cannot be used to further increase the cardinality of the 
matching. Yet, notice that B{[Gq, ... ,G^,I^ ]) may not have a maximum matching of size n; hence, additional 
dedicated sensors indexed by J'” need to be considered, corresponding to the indices of the row vertices that are 
unmatched (i.e., the left-unmatched vertices of B{[Gq, ... ,G'\^,1^']) associated with the MWMM). In summary, 





the indices in J' correspond to the sensors that increase the g-rank and minimize the number of non-accessihle 
state vertices, since they are assigned to non-hottom linked SCCs of T>{Go V ... V Gk), whereas J'” ensures that 
the g-rank is equal to n. Thus, the last step consists in considering an additional subset of dedicated sensors J"’ to 
ensure that there are no non-accessihle state vertices in 'D{Gq V ... V I;^), where J = J' J” U J"'\ in other 
words, J"' consists of indices of the minimum dedicated sensors assigned to a state variable in each non-bottom 
linked SCC that contains non-accessihle state vertices. 

The next result establishes the correctness and analyzes the implementation complexity of Algorithm 

Theorem 4: Algorithm [T] is correct, i.e., it provides a solution to Vi. Furthermore, its computational complexity 
is 0{\C\^), where \C\ denotes the number of column vertices in B{[Gq, ..., G]^, S']). o 

Proof: The correctness of Algorithm [T] follows from noticing that the indices in J' identifies the minimum set 
of dedicated outputs that simultaneously maximizes the increase in the g-rank of [Gj,... ,G^,I^ ] with respect 
to [Gq,...,G^] by \J'\, and the dedicated outputs assigned to state variables in different non-bottom linked 
SCCs. This follows from observing that (by construction) ;B([Gg,..., G^j) incurs in a minimum weight maximum 
matching M with zero weight and size |M|. From Lemma it follows that g-rank([GQ,..., G^j) = |M|. 
Subsequently, a minimum weight maximum matching M' of ;B([Gg,..., G^, 5]) equals \M'\ — |M|; hence, 
increasing by \M'\ — |M| the g-rank of [Gj,..., G^, with respect to [Gj,..., G^j, and contributing to satisfy 
condition (ii) in Corollary [T]- but may not be enough yet to ensure this condition, which is accounted for with the 
set J". In addition, by construction of S it follows that correspond to dedicated outputs that are assigned to 
state variables in different non-bottom linked SCCs; hence, \ J'\ non-bottom linked SCCs have outgoing edges from 
its state variables in different outputs in the system digraph; thus, contributing to satisfy condition (i) in Corollary [T] 
- but may not be enough yet to ensure this condition, which is accounted for with the set J'". Subsequently, the 
number of total number of additional dedicated outputs 1'^" required such that g-rank ([Gj,..., G^, ' j In " ]) — ^ 

is minimized by considering Step 4. Similarly, the number of total number of additional dedicated outputs 1’^'” 
required such that there exist no non-accessihle state vertices in T){Gq V ... V Gj^,is minimized by 
considering Step 4. Notice that 1'^" are not assigned to non-bottom linked SCCs, otherwise they would have been 
considered in 1'^ . Therefore, by setting J = J' ^ J" U J'", as in Step 5, we obtain a solution to V\, since 
Corollary [T] yields. 

The computational complexity follows from noticing that Step 2 can be solved using the Hungarian algorithm that 
finds a MWMM in G(max{|C|, |7^|}^), whereas all other steps have linear complexity; hence. Step 2 dominates the 
final computational complexity, leading to the final complexity of 0(|Cp) obtained, where we note that \C\ > 1^1- 


V. AN ILLUSTRATIVE EXAMPLE 

In what follows, we consider a brain activity dataset |[30ll describing the dynamics of the brain while the individuals 
perform different motor and imagery tasks. The brain activity was measured using the brain-computer interface and 
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consists of 64-channel electroencephalogram (EEG). Each individual performed 14 experimental runs consisting 
of one minute with eyes open, one minute with eyes closed, and three two-minute runs of interacting sessions 
of opening/closing the corresponding left/right fist as a function of where the targets appear on screens. To this 
end, a geodesic sensor net for EEG was considered (illustrated in Eigure [1]), which sensor location schematics and 
enumeration is presented in Eigure 



Fig. 1. 64-channel geodesic sensor net for an EEG. 

Using the wavelet technique described in ll^ . we estimated the parameters aj (j = 1 ,... ,64) in ([T]), corre¬ 
sponding to the fractional order derivatives governing the dynamics of the monitoring neuronal regions and the 
coupling matrix A that encapsulates the spatial fractal connectivity among different hrain regions. To give some 
intuition of how spread the parameters aj corresponding to the fractional order derivatives are, the estimated values 
range between 0.97 and 1.28. Hence, providing evidence of the fractional order dynamics; more precisely, some 
brain regions exhibit pronounced long-range memory dynamics which can be better captured by fractional order 
state space representation. In Eigure we contrast the recorded and simulated data on channel 64 (depicted in blue 
in Eigure [^, during the first two-minutes run of interacting sessions of opening/closing the corresponding left/right 
fist as a function of where the targets appear on screens. 

A common approach in determining the brain connectivity is to consider the correlations between different 
measured signals j^ . In contrast, we advocate for a structural approach of the signal dynamics induced by 
the sparsification of the coupling matrix A. This structural analysis offers a robust mathematical framework for 
capturing both the spatial and temporal interactions among distant brain regions and identifying the minimum 
number of sensors that provide full observable description of brain dynamics. Eollowing the theoretical results 
and analysis outlined in Corollary 2, by resorting to Algorithm 1, we determine the optimal number of sensors 
as a function of the sparsity present in the coupling matrix A for a time horizon equal to the dimension of the 
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Fig. 2. Sensor distribution in the 64-channel geodesic sensor net for an EEG depicted in Eigure The sensor in blue represents the 
sensor whose neuronal activity is simulated using the fractional order system identified, and compared with data recorded in Figure The 
sensors in red represent the minimum number of sensors required to ensure structural observability of a fractional order system identified, 
and considering a sparsity of 80%. 



Fig. 3. In this figure, we show the recorded and simulated data on channel 64 (depicted in blue in Figure]^, where a sample is taken every 
second, during the first two-minutes run of interacting sessions of opening/closing the corresponding left/right fist as a function of where 
the targets appear on screens. 
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Fig. 4. Number of required sensors to ensure structural observability as a function of the sparsity. 

State space. Figure shows that the higher the sparsity of the coupling matrix A of the hrain dynamics is, the 
more sensors are required for ohservahility purposes. An important observation one could make is that for the case 
when the system connectivity is within the range of 20 — 30% (which implies that the sparsity is within 70 — 80% 
interval), which is characteristic to structuraVfunctional hrain connectivity 1(3^ . l(34ll . ll35l . the required number of 
sensors exhibit a first order phase transition. Simply speaking, to monitor the brain activity when the connectivity 
is between 70% and 80% we only need 2 sensors, depicted in red in Figure However, when the connectivity is 
just 10% we are required to use 7 sensors and the required number of sensors spikes fast with increasing sparsity. 

VI. CONCLUSIONS AND FURTHER RESEARCH 

In this paper, motivated by the problem of minimum placement of probes in the brain to obtain dynamic 
observability, we presented methodologies to determine the minimum dedicated sensor placement for discrete time 
fractional-order systems. Eurther, the results presented can be used to determine the minimum dedicated input (i.e., 
inputs that actuate single state variables) placement for aforementioned systems by invoking the duality between 
controllability and observability in the case of discrete-time fractional-order systems. Einally, it would be interesting 
to generalize the present results for a larger class of systems, and validate the proposed methodology in different 
brain regions and using different technologies. 
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